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lOBEiViiaRD 


The  First  Conference  of  Arsenal  Mathematicians  was  held  29  October  195^ 
with  Vfeitertown  Arsenal  acting  as  host.  Colonel  A.  P.  Taber,  CoDonanding 
Officer  at  Watertown  Arsenal,  said  words  of  greeting  to  the  fifty-six  members 
of  the  Conference.  Colonel  P.  K.  Oillon,  Ordnance  Corps,  Commanding  Officer 
of  the  Office  of  Ordnance  Research,  served  as  chaixman  of  the  meeting. 

Colonels  Oillon  and  Taber  pointed  to  specific  ordnance  projects  needing  a 
hi^  level  of  scientific  knowledge,  adding  that  anplled  mathematics  plays  a 
vital  role  in  many  of  these ,  and  advancements  in  this  field  are  needed  if 
hl^  level  improvements  of  ordnance  needs  are  to  continue. 

With  the  one  exception  of  the  Invited  address  by  Professor  F.  J. 

Murray  of  Columbia  University,  the  talks  were  given  by  fhll^tiine  employees 
of  six  different  government  installations.  All  the  papers  covered  points 
of  common  interests  to  those  in  attendance  and  were  of  such  hi^  natheoatioal 
cediber  that  a  general  request  was  made  to  have  the  mazaiscripts  bound  in  a 
single  researdi  report.  In  addition  to  affording  an  opportunity  to  the 
speakers  to  present  their  new  results,  the  conference  permitted  those  in 
one  arsenal  to  see  and  discuss  the  problems  faced  bnr  those  in  another  arsenal. 
Such  Interplay  of  ideas  followed  by  ancxll  rroir  discussions  is  vital  to 
continued  new  attacks  on  the  many  problema  facing  arsenal  mathematicians . 
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of  the  Office  of  Oxdsance  Beaearch 

MCWniS  SBSSION  0900  -  1190 

Introductory  Bemarks 

Stress  Analysis  in  Snail  Arms  Barrels  (30  minutes) 

Br.  A.  B^er,  Springfield  Arsenal 

Ibeory  of  Boise  Modulation  (30  minutes) 

Br.  P.  B.  Karr,  BLamond  Ordnance  lyise  laboratories 

Break 

Influence  of  Elastic  Bsfonnation  on  Accuracy  of  Hl^  BsrfoxaaDce 
Bockets  (30  minutes) 

Br.  H.  J.  Stewart,  Jet  Propulsion  laboratory 

file  Application  of  Electronic  Analog  Cooputers  to  Ptoblass 
of  Tla^  resign  (I9  minutes) 

Mr.  T.  A.  Mood,  Detroit  Arsenal 

Imrii  1150  -  1250 

Ibur  of  Matertown  Arsenal  12^)  -  135^ 

Ammm  session  1330  - 1613 

Problems  in  the  Arolioation  of  least  Squares  to  IVinetioa 
Generator  Design  \y)  minutes) 

Dr.  J.  C.  fii^rt,  Ptankford  Arsenal 

Sboe  T|jrpical  Problems  Confronting  Mathematicians  at  the 
Ballistic  Bsseardi  laboratories  (30  miautes) 

Dr.  S.  Com,  Aberdeen  Proving  Ground 

Break 

Oomputation  Study  of  Differential  Bquations  (30  minutes) 
Professor  f.  J.  Murray,  Oolonbia  Ibiveraity 
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Vbtertown  Arsenal  29  Octoter  19^ 

Colonel  P.  N.  Oillon,  Ghalxmazi 

siRBiss  nr  sull  aims  barrels 

I|jr  A.  Banmer 

Increased  firing  rates  and  the  ever  Increasing  severity  of  the 
firing  schedules  induce  conditions  in  'barrels  of  snail  arms  weapons, 
from  Cali'ber  .30  to  3Qb>>&i  ^dilch  require  different  design  and  analytical 
approaches  from  those  used  in  the  past.  Vital  points  of  new  harrel 
constructions  are  descrihed  and  discussed  and  'brief  outlines  of  the 
processes  %hich  occur  within  lined  barrels  are  given  in  the  first  part 
of  the  paper.  The  barrel  stress  analysis,  vhich  is  the  subject  of  the 
second  part,  considers  the  themal,  shrink-fit  and  rifling  stresses  as 
well  as  the  pressure  stresses,  and  the  graphical  presentation  of  the 
total  equivalent  stresses  facilitates  the  task  of  establishing  optiaom 
design  parameters  for  ll^twei^t  accurate  betrrels  of  extenAed  life. 


The  main  portion  of  the  address  by  Dr.  A.  Hemner  can  be  found  in 
the  security  information  paper  entitled  "Barrel  Research,  Development 
and  Stress  Analysis  at  Sjpringfield  Armory."  Cbpies  of  this  paper  may  'be 
olttained  from  fljpringlleld  Armory  by  authorised  persons. 
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A  NOTE  ON  ZERO  CBosscTos  iun)  AvimoE  mpiNcr 

F.  H.  Karr 


1.  Introduction 

The  function  alnCRffg^t-*-  ^T)  with  the  constant  freqaanegr  io»  the 
constant  angle  IS.  t  and  t  denoting  ^  tine  varlablef  bepenea  aero 
2g0  times  per  unit  time.  A  natural  qjueatioa  that  arlaaa  las  if  in 
the  above  expression  the  phase  funotion  2ag^t  ia  raplaead  bgr  a 

random  function  of  time,  at  what  average  rata  will  aaro^oreaeinga 
occur  and  to  what  extent  can  thLa  rata  be  aaaoaiatad  with  an 
''average  frequency"  of  the  random  funotjont  Ihia  qiuaatlom  la 
discussed  In  the  following. 


2.  Development 

We  consider  first  the  funotion  8ln(N^t‘)v  if )  where  ^  la  a  oonatant 
as  above  and  N(t)  Is  a  stationary  random  ("noise")  funotion  of  t|  In 
the  present  paper  we  assume  that  N(t)  has  a  Gaussian  amplitude  dlatrl* 
butlon.  To  find  the  average  number,  n,  of  saroa  of  aiB\ll(t)a  if  )  m 
first  find  the  average  number  of  tinas  par  unit  time  that  the  arpnant 
N(t)+  0  crosses  the  values  0,  ^2*,  ato}  we  then  add  all  thaaa 

results  together,  giving  an  infinite  series,  whose  sum,  it  ia  aaaa,  is 
the  desired  result,  n. 

(1) 

The  fundamental  oaper  of  Rice  treats  the  problem  of  the  rate  of 
zero-crossings  of  N(0.  To  make  the  caloulations  listed  in  the  praeedimg 
paragraph,  we  should  like  a  formula  for  the  ejqwoted  nuidwr  of  *A«ereaalngi*, 
that  Is,  the  expected  number  of  oroaslngs  per  unit  time  of  the  aiUtrarjr 
value  "A"  by  the  noise  function  N(t)«  It  turns  out  that  auoh  a  fesamla 
can  be  obtained  by  a  simple  extension  of  Rioe'a  darimtion  for  the  aere* 
crossings.  The  result  is 


r  ■ 


wh(o) 


1/2 

-1 


-a2 


mo) 


(1) 


where 


r  Is  the  expected  number  per  unit  time  of  oroaainga  of  the  ealma 
"A"  by  N(t), 

R(t)  is  the  oorrelatlon  funotion  of  N(t), 

E(t)  Is  the  second  darlvatlva  of  R(t)  with  respeot  to  t, 
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and  R(o)  and  R(o)  are  tha  values  of  R(t)  and  R(t}  at  t  >  0.  As  is 
well  known«  R(o}  is  the  mean  Intensity  of  N(t),  or  M  • 

Now  N(t)4]r  has  the  value  of  A  vhen  N(t)  has  the  value  of  A» 
so  to  obtain  n  we  need  to  give  A  in  eg*  (1)  the  values  Idr  -V  for  all 
integral  values  of  k,  and  add  the  resulting  terms.  We  have  then 


n  -  c->:L5i£L3 
n^R(o) 


1/2 


exp  • 


(to 


(2) 


Ic  «  •  so 


2R(o) 


Eq.  (2)  is  a  fomula^  to  be  discussed  in  more  detail  below,  giving  n 
in  terms  of  ,  R(o)  and  R(o).  Before  considering  the  dependence  upon 
these  factors  however,  we  tdsh  to  transform  (2)  into  a  fora  which  will 
exhibit  specifically  the  ^average  frequency**  of  N(t)* 


We  havf,  the  Wlener-Klnohine  ttaeorea 
R(t)  -  y  S„(f ) 


(Z) 


cos  2w  ft  df 


(3) 


where  Sjj  is  the  spectral  Intensity  or  *power  spectrum**  of  N(t)  and  f 
denotes”f reque  ney • 


Then 


-  R(t)  - 

j  Uw2f2sj,(f )  cos  2Wft  df* 

U) 

Also 

R(o)  ■ 

}  S,  (f)  df 

(5) 

and 

-  R(o)  - 

0O 

)  l4r^S„  (f)  df* 

S  A  " 

(6) 
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Now  we  may  Introduce  the  **lnstantaneous  frequency*  of  N(t)  bj 
assuming  it  to  be  expressible  In  the  fom 

t. 

N(t)  -  2i^  g(t)<tt  (7) 

£4*  (7)  expresses  the  relation  betf#een  N(t)  and  g(t)  aa  that 
between  phase  angle  and  Instantaneous  frequency*  Since  the  tranefer 
function  of  a  perfect  Integrator  le  ,  we  have 

2«jf 


SK(f ) - -  S,(f) - ^  (d) 

hm  i  t 

where  Sg(f)  is  the  spectral  intensity  of  g(t)* 

Therefore 

Oat 

-  a(o)  .  C  to'  a,(f)  it.  (») 

Now  it  is  a  well  known  property  of  statioiary  randoa  prooeese* 
such  as  g(t)  that 

T  •• 

11«  — i-  (  g^t)  dt  -  C  S  (f )  df  -  g*  (ID) 

i  i  ^ 

where  g^  denotes  the  time  average  (or  stochaetic  average)  of  g^« 

Eq*  (9)  therefore  has  the  interpretation  that  if  a  noise  function 
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N(t)  may  be  expressed  as  the  Integral  of  an  Instantaneoiu  frequency, 
and  If  R(t)  is  the  correlation  function  of  N(t),  then  •  R(o)  Is 
times  the  mean  square  frequency  of  N(t)« 

At  this  point  vie  introduce  the  average  value  of  the  absolute  value 

Isl*  (The  average  g  is  zero*)  If  g(t)  has  a  Qaussian  amplitude  dis¬ 
tribution,  vihieh  ue  assume,  and  vdiich  is  consistent  with  the  pravious 
assumption  that  N(t)  is  Qaussianly  distributed,  then  it  is  easily  found 
that 

« 

Introducing  the  last  tvfo  results  ve  may  now  rewlte  (2)  in  the 

form 


An  applloation  of  Poisson's  Formuli^^o  (13)  gives  the  alternate 


The  task  of  numerical  computation  is  greatly  facilitated  by  these 
two  alternate  forms  since  (13)  converges  rapidly  when  R(o)  is  small, 
while  (111)  is  rapidly  convergent  vdien  R(o)  is  large*  Thus  the  two 
alternate  forms  oomplement  each  other  nicely* 
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From  (lU)  we  the  intereating  result  for  large  R(o)^ 
and  Independently  of  o  >  n  la  asyopto^eally  2|g|y  that  l8«  twloe 
the  average  frequency  (the  oognonen  (^average  freqiosnoy*  applied  to 
Tgl  lapli^*  of  oourae  that  negative  and  positive  frequenoles  are  re¬ 
garded  as  on  an  equal  basis)*  For  small  R(o)  n  Is  aa3mq>totloally 


In  Flg*^  there  Is  plotted  the  ftinotion  Til  for  various  values 
of  R(o)  and  (T  •  It  will  be  seen  that  n  can  be  greater  than  2|g|f  as 
well  as  Issa.  The  Influence  of  If  becomes  rapidly  dissipated  as  R(o) 

Increases*  In  fact  Fig,  1  shows  that  even  for  R^  ■  — iSL.— n  Is  already 
essentially  equal  to  2TgT«  ^ 


One  may  ask«  what  is  tbs  value  of  n  when  averaged  over  ^  2  The 
font  (lU)  Is  oonvenlent  for  examining  this  question.  Ohs  sees  that  all 
the  oos  2yn  terms  drop  out  upon  averaging  and  we  are  left  with 


(15) 


Thus  the  average  value  of  n  is  the  same  as  the  asymptotic  value  for 
high  mean  Intensity  R(o),  namely  twloe  the  average  frequency*  Thus  for 
either  of  these  conditions  the  oonoept  <*average  instantaneous  freqyienoy* 
Is  useful  and  valid,  since  the  expected,  number  of  sero  crossings  Is  twice 
the  average  frequency*  This  is  an  intuitively  satisfying  generallaation 
of  the  corresponding  result  for  sin(2w^  g^dt  *}f]i  that  is  under  the 


conditions  indicated.  If  g^  In  the  preceding  ejqpressi^injM  replaced  bv 
g(t)  then  the  result  n  ■  2gQ  Is  generalised  to  n  ■  2|g)  or  n^^^  ■  2rg|  * 


3*  Suitinary 

An  e:iq;)resslon  Is  derived  for  the  expected  rate  n  of  sero  orossli^s 
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of  the  function  8tntN(t)  *  ^)t  where  N(t)  is  a  stationary,  Qausslanly 
distributed,  random  or  ’*nolse<*  function  with  zero  mean.  The  e:qpresslon 
for  n  Involves  n(o),  the  second  derivative  with  respect  to  t  of  the 
correlation  function  R(t)  of  N(t)  at  t  ■  0,  R(o),  and  .  When  N(t) 

Is  expressed  as  the  Integral  of  an  Inst^taneous  frequency  function 
g(t).  In  the  form  ll(t)  ■  2*^  g(t)dt,  R(o)  may  be  expressed  In  terms  of 

nr,  the  average  value  of  |g|.  ^le  resulting  e^qpresslon  for  n  Is  pro¬ 
portional  to  |g| ,  the  proportionality  factor  being  a  function  of  ^ and 
R(o),  the  mean  intensity  of  N(t).  When  R(o)  Is  large,  this  factor  Is, 
closely  eq^al  to  2  for  any  value  of  ,  while  for  any  value  of  R(o), 
the  value  of  this  proportionality  factor  averaged  over  all  values  of  0 
Is  exactly  2«  These  msults  are  In  pleasing  analogy  with  the  case  In 
whldi  g(t)  Is  replaced  by  the  constant  g^.  In  tddch  case  n  Is  2gQ.  Thus 
the  concept  of  ^average  Instantaneous  frequency**  is  a  useful  one  In  the 
sense  that  the  ratio  of  2  that  holds  between  n  and  the  fixed  frequency 
g.  also  holds  between  n  and  the  averam  frequency  ig|  If  either  of  th^ 
above  conditions  is  considered.  If  Rto)  Is  small  and  averaging  over  0 
Is  not  allowed,  the  sltuatiim  Is  less  simple.  Details  are  given  In 
formulas  and  In  Fig.  1. 

U»  Acknowledgements 

This  work  was  stimulated  by  discussions  with  Hr.  B.  N.  Horton  tAio 
has  been  interested  for  some  tlms  In  the  concept  of  **average  frequencjp*. 

Most  of  the  numerical  work  associated  with  Fig.  1  was  performed  by 
Mr.  C.  Morrison. 


S.  0.  Rice,  Bell  Telephone  System  Monograph  1569,  p.  57* 

N.  Wiener,  Aoto  Math.  $$,  117  (1930). 

Morse  and  FMhbaoh  "Nsthods  of  Theoretical  Fhyslcs*  p.  U67. 


s 


OonfersiiM  of  Araenol  Matheoatieiaat 


Conference  of  Arsenal  Matbanaticlans 


9 


immiENCE  OF  msnc  iefoimation 

ON  ACCDSACT  OF  HIGE  PBRFOMAIICE  BOGKETS 
H.  J.  Stewart 


The  paper  presented  hy  Dr.  H.  J.  Siteiart  contained  some  material  abstracted 
from  reports  having  a  Confidential  classification.  Those  interested  in  porsoing 
further  the  subject  matter  he  discussed  are  referred  to  the  following  articles 

1.  JFL  Publication  Ho.  24.  ^  J.  Z.  Froehllch  (1933)  • 

2.  JPL  Progress  Beport  Ho.  20-290.  by  I.  Statler  (In  publication). 

3.  JEL  Progress  Beport  Ho.  20-251,  b7  J*  loarell  (In  publication). 
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m  USE  OF  mi£>a  CQtroiXRS 
IN  IBE  IZSSKBf  OF  MILXTABI  HlBICiaB 
ly  fhomae  Wood 

!Q)e  development  of  militaxy  tanks,  like  other  auto¬ 
motive  vehicles,  depends  not  only  on  mathematical  ccnputa- 
tion  hut  also  on  a  procedure  Ixised  on  construction  of  models 
and  extensive  teste.  It  has  not  been  unusual  to  construct 
complete  vehicles,  subject  these  to  severe  laboxatoiy  and  ardu> 
ous  road  tests,  only  to  discover  that  the  new  design  may  not 
have  contributed  materially  to  the  desired  perfozmanee  character¬ 
istics  of  the  vehicle.  Add  to  this  a  shortage  of  trained 
personnel  capable  of  conducting  a  mathematical  analysis  and  it 
can  he  readily  understood  costs  of  sagiBssrlag  dsvalni— it 
nay  grow  out  of  proportion. 

Blanks  to  electronics,  relatively  low  cost,  hi^  qpsed 
analog  computers  have  became  available.  Bigineers  at  the  Det¬ 
roit  Arsenal  recognized  the  possibilities  of  this  instrument 
and  now  utilize  the  computer  in  the  eetrly  stages  of  the  develop¬ 
ment  of  a  complete  vehicle  or  its  complex  components.  The 
advantages  %diich  result  from  this  rapid  analytical  investiga¬ 
tion  over  the  time  consvming,  tedious  and  costly  manoal  com¬ 
putations  are  numerous,  but  more  important  is  the  fact  that 
new  ideas  and  desigi  principles  can  be  investigated  almost  as 
xs^idly  as  they  are  conceived. 
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The  computer  makes  it  possible  to  economically  detexmlne 
the  feasibility  of  a  design  vdthout  constructing  expensive 
models  and  conducting  extensive  laboratory  teste  vftilch  are 
often  destructive.  Moreover,  it  helps  the  engineer  to  de¬ 
termine  the  design  principles  %iiich  offer  the  most  promlsii^ 
avenues  of  investigation.  In  cases  vhere  uncertainty  exists 
concerning  the  values  of  parameters,  the  computer  aids  in  de¬ 
fining  the  limits  within  \dildi  the  parameters  can  vary  to 
produce  satisfactory  results.  This  feature  is  not  only  advan¬ 
tageous  in  desi^  and  development,  but  is  also  of  extreme 
value  in  establishing  tolerances  on  vehicle  parts  idiich  have 
been  approved  for  production.  This  prevents  the  costly  manu¬ 
facture  of  precision  parts  for  a  system  in  vhich  precision  is 
not  necessary  for  performance. 

An  analog  computer,  as  its  name  implies,  is  tan  electronic 
instrument  capable  of  makiiig  hi^  qpeed  mathematical  calco- 
latlons,  once  the  x^slcal  b^vior  of  a  device  has  been 
described  and  formulated  into  mathematical  statements.  This 
can  be  done  because  the  behavior  of  most  physical  systems  can 
usually  be  expressed  in  terms  of  mathematical  equations  and 
grc^s.  Mhen  the  mathematical  or  ^if^ical  statements  become 
available  and  the  various  elements  (multipliers,  integrators, 
etc.)  of  the  computer  are  connected  in  a  manner  analogous  to 
the  equations,  the  computer  will  Tjerform  the  appropriate  naths- 
natioal  operations. 
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!Ihe  more  frequent  f\jnclaaental  operations  perfomsd  Igr 
analog  computers  are: 

1.  Multiplication  ly  a  constant. 

2.  Integration  of  qua.  titles  with  rest:>ect  to 
time. 

3.  Multiplication  of  varlahles. 

Cie  operational  asipllfier  is  the  lasle  huilAlng  hloek 
of  the  computer.  The  mrious  operations  t^ileh  it  performs 
depends  on  the  Impedances  that  are  placed  in  the  fOedhaek  and 
input  circuits  of  the  amplifier.  All  the  phgrsioal  panmetere 
are  expressed  in  terms  of  roltages. 

TO  Illustrate  the  method  hy  which  the  operational 
amplifier  of  the  computer  functions  as  a  multiplier,  let  um 
assume  that  (Heart  X)  the  pdn  0  of  the  amplifier  appromehma 
Uiflnitgr. 


A 

e- 


AivuniR 


nOURI  I  -  MULTWLICATIOM  BY  A  COMTAMT 
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To  get  a  finite  voltage  of  the  amplifier  at  B  vie  mast  not 
have  any  voltage  at  point  0;  hence  Eg  vdll  he  zero.  In  a 
practical  case  vie  may  use  a  gain  in  the  order  of  200,000 
vdiich  makes  Eg  not  zero  hut  of  such  magnitude  that  it  can  he 
Defected.  By  placing  resistors  R2  and  in  the  feedback 
and  input  circuits  respectively,  the  amplifier  is  said  to  he 
"patched"  for  multiplication.  As  an  example,  let  us  set  the 
operational  amplifier  to  multiply  ty  tvo  "ty  placing  a  one 
megohn  resistor  fji  the  feedhstck,  and  a  .3  megohm  resistor  in 
the  input.  Ihe  ope  ration  is  performed  because  aic^  voltage  ^ 
inpressed  at  A  across  che  input  resistor  B]^  causes  a  current 
II  to  flovi  in  the  irput  circuit.  A  ftodamental  electrical 
theorem  (Xlrchhoff's  lav)  states  that  the  algebraic  sum  of  the 
currents  that  meet  at  any  jiuK:tion  is  zero.  Therefore 

^2  ^  Ig  must  equal  to  zero.  Ig  is  the  current  drawn  ’qy 
the  grid  of  the  operationeu  amplifier  and  is  on  the  order  of 
10"^  ampere-,  therefore,  i:^  cem  he  neglected.  Thus  I2  must 
he  equal  in  magnitude  and  opposite  in  sign  to  to  make  the 
son  equal  to  zero.  This  produces  a  sign  change  across  tbs 
aopllfier  and  will  induce  a  voltage  which  is  of  opposite 
sign  aad  twice  This  is  tree  because  I2  is  equal  to 
and  I2  ie  twice  lEt ,  Hence,  is  tvdee  To  multiply  1y 
asy  constant,  it  is  necessaxy  to  select  the  appropriate 
(Bg  is  usually  kept  at  1  megohm) . 
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Tigare  2  illustrates  the  method  of  "patching"  the  operatlonftl 
amplifier  for  use  at  an  integrator. 


rwiffit  1  .  MTiaunoN 


"Bf  jdacing  a  1  neigohn  reeletor  in  the  input  eizcuit  (1^) 
and  a  one  aioroflaxad  c^paeator  (0^)  in  the  flMdlack  oireuit, 
the  opemtional  anplifler  will  genezate  the  tine  integzml  of 
an  input  ^Itage  Again  using  SLrdioff'  e  lav  we  eee 
that  the  woltage  at  A.  called  produeee  a  current  to 
flow  throu^  resistor  In  order  to  satisfy  the  law  that 
the  algehraic  sun  of  the  currents  %Aiich  neet  at  angr  Junction 
is  equal  to  serot  X2  ^  equal  in  magnitude  and  opposite 
in  sign  to  I^.  To  produce  this  current  Xg  the  voltage 
across  O2  must  he  increased  at  a  omstant  rate;  hence  1^  is 
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the  lnto(5ral  of  E^, 


nsure  3  3ho«  that  th.  alop.  of  1.  depoodent  upon  tho  ratio 
Of  Oj  and  Hj.  Tharafore,  »«  ana  both  Intaerata  and  nultlply 
by  a  constant  with  this  arruncenent. 

amtlplloutlon  Of  t«o  Turlablaa  oan  b.  porfomad  b,  nawia 
Of  a  sonro  system. 
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I, 


noMt  4  .  iMthwue*imi  cr  nM  vmimim 


la  Ti^gaan  4  the  aeviatle  prodooee  aa  oatpat  deeft  aatlaa 
whoee  aagilar  poeltion  -mrlee  id.th  respect  to  time  escaetlj 
as  Ths  oa^t  shaft  controls  the  position  of  tbs  an 
of  a  potentloaeter  aoxoss  «hidi  the  variable  Bg  is  applied. 
Tbs  Tolta^  on  the  aim  is  a  ftaetioBal  part  of  ths  applied 
eolta^  X2t  ^  ■  ^^2* 

Althoo^  ths  possibility  of  using  an  analog  eoaputer 
for  design  and  dsTslopaent  of  military  vehicles  was  appeirsnti 
there  was  no  way  of  dstamining  the  aocoracy  inherent  in 
this  peurticolar  ^rpa  of  simulation.  Ibr  this  reason,  a 
research  progxam  was  instituted  at  tho  Willow  Ron  Beseardi 
Center,  under  contract  with  ths  Detroit  Arsenal,  lbs 
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purpose  of  the  research  vas  to  correlate  the  reQ)onse  of  an 
actual  tank  suspension  system  with  that  of  a  simulated  tank. 

AnM47  tank,  welding  approximately  4S  tons,  was  select¬ 
ed  for  experiments  in  order  to  correlate  the  simulated  con¬ 
ditions  with  the  experimental  data.  The  features  of  the  taiic 
are  shown  in  Figure  3* 


fwnt 


The  bull  ie  suq[)ported  ly  twelve  road  viwels,  six  on  es^h  sidbs, 
Qisse  ilisels  are  attadisd  to  the  hull  throo^  road  aras. 


IS 
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shock  absorhers,  and  torsion  bars  which  in  effect  act  as  springs. 
Volute  springs,  or  bump  stops,  as  they  are  commonly  called,  are 
installed  to  limit  the  deflection  of  the  wheel  assembly  vdth  respect 
to  the  hull .  Attached  to  road  vdieels  1,  2,  ^  and  6  are  shock 
absorbers  vhose  function  is  to  damp  out  motion  of  the  hull.  Tb 
increase  the  mobility  of  the  tank  over  rou^  terrain,  a  track  is 
provided,  on  eadi  side  of  the  tank. 

(Bib  field  tests  consisted  of  running  the  tank  at  varions 
speeds  (2  to  2D  mph)  over  two  types  of  road  obstacles.  One 
obstacle  M&e  a  12**  and  12"  timber,  and  the  other  an  Inserted 
7  ramp  10'  loi^  and  8"  hl^,  as  ^%fn  in  figure  6. 


V  .iil... 


fii.  I: 


-.1 


NOTES. 

1.  PoiHiv*  diractlM  of  fofco  and  dliptacomant  li  auumod  upword.  PotMvo  dirocHon 
of  momonf  and  angla  It  auumod  doctnuita. 

3.  Solnranca  podtion  of  hull  and  ooch  whool  It  vortleai  potWan  of  obeh  Mom  If  tank 
woro  rotlint  on  lowol  ground  and  hod  no  wolght. 

Raforonco 

ghiililrf 

ofNali 
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Seta  pertinent  to  the  motion  of  the  idieele  and  boll  viere 
recorded  on  motion  pictures  trca  which  the  pitch  an^e  and 
vertical  dlQ)lacement  of  the  hull  and  each  of  the  road  whaels 
could  he  measured. 

Sven  a  cureozy  analysis  of  the  functioning  of  this 
suspension  system  makes  it  obvious  that  to  take  into  account 
all  effects  would  require  a  computer  setvqp  of  formidable 
else  with  increased  time  and  cost,  fo  keep  the  simulation 
within  reasonable  limits,  it  was  therefore  desirable  to 
look  for  simplifying  assumptions  whidi  would  not  seriously 
affect  the  accuracy  of  the  results.  Ito  accomplish  this,  the 
tank  was  considered  as  a  solid  body  in  space  which  has  two  de¬ 
grees  of  freedcxn.  The  two  degrees  of  freedom  which  were  ooa- 
sidarsd  in  the  slanlatiea  of  tbs  NUJ  tank  for  the  analOf 
cooiputer  were: 

1.  7ertical  motion  of  the  center  of  gravity 
(usually  called  bounce). 

2.  Botatlonal  motion  (usually  called  pitch)  about 
a  transverse  axis  running  horisont^ly  throu^ 
the  center  of  gravity  of  the  tank. 

In  the  simulation  (Figure  7)  >  the  hull  of  the  tank 
considered  to  be  a  rigid  mass  coupled  to  six  masses  (wheels) 

'  ^*V  six  sprii^  (torsion  bars  plus  volute  springs)  and  four 
dampers  (shock  absorbers).  These  six  masses  were,  in  turn, 
si^ported  above  the  road  by  six  springs  and  dampers  represent¬ 
ing  the  cpring  and  dasiping  effect  of  the  rubber  tires  and  track. 
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Only  six  vftieele  and  half  the  mass  of  the  tank  were  considered 
since  it  vas  assumed  that  the  tank  is  symmetrical  ahoat  a 
longitudinal  vertical  plane. 

In  general,  it  is  convenient  to  start  with  an  over¬ 
simplified  simulation  and  make  further  refinements  with  their 
attendant  complications  as  the  need  hecomes  evident.  The 
major  assumptions  which  were  finally  found  to  he  satisfactory 
were  not  all  detezmined  in  advance,  hut  grew  out  of  the 
correlation  work  itself.  Careful  observation  of  the  test 
results  gave  clues  to  the  factors  which  diould  not  he 
neglected  for  realistic  simulation.  For  examule,  it  was  con¬ 
venient  in  the  initial  comnuter  studies  to  assume  that  the 
vteels  did  not  leave  the  ground.  However,  observations  from 
the  motion  picture  records  clearly  estahlidiod  the  inportaooe 
of  simulating  this  effect. 

The  assumptions  concluded  and  used  are  as  follows: 

(l) .  The  track  has  only  minor  effects  on  the  hull 
response  and  if  necessary  may  he  neglected  in  the 
simulation.  The  only  effect  of  the  track  included 
in  the  simulation  was  the  modification  vdiich  results 
from  the  action  of  the  track  as  it  affected  the 
vertical  motion  of  the  road  wheels  when  they  moved 
onto  and  off  the  timber. 

(2) .  The  spring  effects  of  the  torsion  bars  and 
the  hump  stops  act  in  a  vertical  direction  only  and 
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are  represented  'by  a  corre  of  force  versas  deflection 
Bide  up  of  a  series  of  stral^t-line  se^Mats  as 
lAaown  in  TLgare  S. 

(3)  •  The  shock  ahsorher  chaxsicterl sties  can  also 
'be  represented  'by  vertical  cooponents  only,  using 
curves  made  of  a  series  of  stxal^t-line  sagaents 
as  also  shown  in  Tigure  8. 

(4)  .  The  horisontal  fbrees  on  the  tank  can  he  neg¬ 
lected. 

Ihe  anolog  computer  detexmined  the  motion  of  the  hull  in 
pitch,  and  the  vertical  deflection ,  as  fhnetions  of  time  by 
making  uee  of  equations  expressing  the  foUotring  aatheaatioal 
xalatloae. 

(1) .  Hhe  son  of  the  vertical  forces  on  the  hull  is 
equal  to  the  mass  of  the  hull  times  its  vertical 
aocelexatlon.  Qie  vertioal  forces  are  ttiose  due 

to  the  torsion  hare,  the  volute  springe,  the  diode 
ahsorhere,  and  gravity. 

(2)  .  The  sum  of  the  moments  about  the  pitch  cuds 

thxou^  the  center  of  gravity  is  equal  to  the  moment 
of  inertia  of  the  hull  with  respect  to  its  center  of 
gravity  times  tbs  acceleration  of  the  hull. 

(3)  •  Ihe  deflection  of  each  road  dieel  with  respeet 
to  the  hull  is  a  Ihnction  of  (a)  the  vertical  positloa 
of  the  road  whsal,  (h)  the  vertical  position  of  the 


Vartleal  motion  o(  tho  Imll.  7.  Vertical  motion  of  wheel  6. 
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FIf.  11  COMPARISON  OF  AAEASUREO  AND  SIMUUTEO  VALUIS  OF  PITCH 
AND  VERTICAl  MSPUCEAAENT  OF  HUU  UNOIt  STATIC  CONOITIONI 
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center  of  gravity  of  the  hall,  (c)  the  pitch  an^e 
of  the  hull,  and  (d)  the  knowfi  geometry  of  the 
soGpenelon. 

(4).  The  sonmation  jf  vertical  forces  on  each  road 
vftieel  Is  equal  to  the  mass  of  the  vftieel  times  its 
vertical  acceleration.  (Che  vertical  forces  on  each 
wheel  are  those  due  to  the  torsion  har,  the  volute 
spring,  the  shock  ahsorher  and  the  reaction  due  to 
the  spring  and  the  danplng  properties  of  the  mhhar 
in  the  track  end  tire. 

The  equations  resulting  from  the  above  matfaematioal 
latlons  are  diovtn  in  figure  9»  ^  figure  10  Aam  the  sche 
■mtic  diagram  of  the  vehicle  simulation. 

▲s  a  means  of  oomparing  the  real  tcuk  with  the  steHatit 
tank,  under  the  least  complicatea  conditions,  cheoks  were 
made  of  hull  pitch  and  vertical  deflection  under  etatic  oou> 
ditlons.  These  data  were  obtained  by  riding  the  tank  unto 
an  obetacle,  allowing  it  to  come  to  rest,  and  then  meatsnring 
pertinent  dietances,  from  which  vertical  dieplacemant  and 
pitch  of  the  hull  could  be  detezminod.  Oomparing  these 
maasured  quantities  with  the  simulated  values  provides  a 
dwck  on  the  validity  of  the  oooputing  technique  for  the  steady- 
state  conditions.  This  comparison  is  shown  for  the  ease  of 
the  l^inoh  square  obstacle  in  Tigure  11,  where  pitch  and  verti¬ 
cal  displacement  are  plotted  Ibr  dlffszwnt  horlaontal  looatloBs 
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of  the  tank  vdth  respect  to  the  hump.  It  con  he  seen  that 
the  correlation  between  real  and  simulated  values  is,  for  the 
most  part,  within  the  limits  of  the  measurement  errors. 

An  e^Eonple  of  the  correlation  of  the  dynamic  response 
of  the  tank  for  runs  over  the  inverted  7  raa^  at  20  fps  is 
Shown  in  Figure  12.  Hie  correlation  of  the  response  of  the 
tank  for  runs  over  the  12"  square  timber  is  also  diown  in 
Figure  12. 

In  addition  to  the  work  on  tank  suspension  systans , 
oonsiderahle  effort  has  been  esqiended  on  the  simulation  of 
tank  gun  control  systems  and  components.  Special  attention 
has  been  given  to  the  problem  of  stabilizing  the  tank  gun  in 
space.  It  has  been  possible  to  predict  with  reasonable  accuracy 
the  performance  of  a  given  gun  stabilization  system  vhen 
coupled  to  the  aforementioned  suspension  simulations.  Hie 
addition  of  a  pictorial  presentation  of  target  and  si^t-fleld, 
along  vdth  the  necessary  controllers,  permits  the  human  gunner 
to  be  included  in  the  overall  simulation.  A  simulation  sudi 
as  this  allows  the  vehicle  to  be  evaluated  as  a  complete, 
dynamic  fitting  unit.  Such  an  evaluation  is  of  much  greater 
value  than  attempting  to  evaluate  the  individuetl  components 
alone. 

Another  problem  in  vhlch  the  analog  computer  has  been 
usefbl  is  in  the  suspension  of  hl^-inertia  guns  and  turrets 
on  moving  hulls.  In  cases  where  gun  or  turret  inertia  is 
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REAL  vs.  SIMUUTEO  RESPONSE  OP  M47  (•  Inch 
Invnrtad  •  V  Ramp,  Tank  Sp«*d  of  20  fps.) 


—  Slmwlornd  Dwa 

Fig.  12  REAL  vs.  SIMULATED  RESPONSE  OP  M47  (12  ktch  Squorp, 

Tank  Sp—d  of  11  fpc.) 
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extremely  it  'becomes  difficult  to  "hold"  the  gun  or 

turret  relative  to  the  hull  vfeile  moving  over  rou^  terrain. 
lEhie  problem  has  been  overcome  by  the  ue»  of  resilient  mem¬ 
bers  ‘between  hull  and  gun.  The  design  of  these  resilient 
members  is  closely  tied  in  with  the  vehicle  suspension  and 
gun  control  r^stem.  Simulation  of  the  suspension,  gun,  and 
resilient  mount  has  allowed  the  detezminatlon  of  design  pare^ 
meters  vMch  would  otherwise  have  had  to  'be  detexmlned  fixsm 
vehicular  tests. 

Ixperimental  work  is  also  in  progress  to  provide  the 
hunan  gunner  a  movable  seat  arrangement  which  will  acc^t  in¬ 
puts  frm  the  simulation  and  impart  the  vehicle  motion  to 
the  gunner.  It  is  recognized  that  the  effectiveness  of  the 
gunner  (and  other  crewmembers)  is  affected  by  vehicle  motion. 
IBie  ability  to  study  these  effects  under  controlled  conditions 
will  ‘be  of  great  value. 

^so  planned  is  a  function  generator  to  impart  a  realistic 
terrain  function  to  the  simulation,  ^topographical  maps  of 
various  terrains  in  different  parts  of  the  world  where  tank 
warliBre  is  contemplated  would  be  reproduced  and  studied. 

A  word  of  caution  may  be  appropriate.  Ihe  analog  com¬ 
puter  is  an  InstiTment  capable  of  solving  numerous  simultan¬ 
eous  equations  which  result  from  effectively  isolating  each 
component  system,  and  usually,  the  equations  for  such  component 
systems  are  of  the  second  order  and  the  number  of  these 
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equations  is  of  no  consequence  "because  the  computer  will  eliminate 
the  time  necesssary  to  solve  them.  Hb^irever,  the  application  of  an 
electronic  computer  to  the  solution  of  an  erigineering  prohiem  depenls 
on  the  availa'bllity  of  numerical  data.  Hence,  effectiveness  of  the 
canputer  approach  will  depend  on  the  completeness  and  validity  of 
these  data.  This  is  not  as  disadvantageous  as  it  may  sean.  Let  us 
assume  that  a  radically  new  eystem  is  to  "be  analyzed,  that  the  Infozna- 
tion  may  "be  incomplete,  and  requires  extensive  extrapolation.  In  this 
case  the  computer  will  provide  qualitative  results  as  to  the  amount  and 
manner  in  which  each  design  variable  affects  the  system's  ^rfoznance. 

The  simulation  will  indicate  the  design  of  the  system  whidi  has  optimam 
characteristics,  "but  will  net  provide  an  absolute  measure  of  perfoxnance. 

Ihe  foregoing  research  j)rogiBm  verified  that  the  results  obtained 
"by  the  campiter  were  valid  and  vdthin  limits  for  the  design  and  develop- 
ment  of  probleme  Involved  in  various  types  of  wheeled  and  tracked  v^cles. 

Wjlle  it  is  true  that  computer  tedmlquss  will  save  an  exsozmous 
amount  of  time  and  money  x^ch^^zmally  would  "be  qpent  in  construction 
cost  and  lahoratozy  teste,  mucdi  more  could  be  accomplished  if  more 
engineers  were  adept  in  corrslatlr,g  the  physical  reaction  of  structures 
in  tezms  of  gzaphlcal  data.  This  difficulty  was  encountered  and  vas 
solved  "ty  developing  a  "simuscope".  This  is  a  visual  aid  vhich  pre¬ 
sents  a  reproduction  to  scale  of  the  motion  of  the  hull  and  wheels 
as  the  veliicle  goes  over  an  obstacle  at  different  speeds. 

In  view  of  the  claims  vMdi  have  "been  made,  one  venders 
computers  are  not  used  more  A:^uently.  The  reasons  are: 
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1.  A  lade  of  loaowled^  of  just  viiat  computers  can  do  to  solve 
proldems  in  a  particular  engineering  field.  The  automatic 
computing  field  is  a  relatively  new  one  and  the  "word"  has 
not  been  passed  on  to  everyone. 

2.  A  mistaken  idea  as  to  the  cost  of  computing  equipment. 
Ihgineers  vho  have  not  had  time  to  become  familiar  with 
computers  feel  t})s.t  it  takes  ^  million  dollar  machine" 
to  solve  problans.  This  is  not  the  case  at  all.  If  one 
were  to  hazaxd  a  guess  on  the  cost  of  a  machine  to  solve 
a  particular  problem  it  vnuld  probably  go  something  liloe 
this:  vhere  computers  are  applicable,  the  cost  of  the 
commuting  equipment  would  be  much  less  than  one  costly 
"guestimate"  vhere  the  entire  design  has  to  be  scrapped 
after  the  pilot  model  is  mamlhctured,  and  in  many  cases 
the  cost  is  less  than  dhanglng  major  components  in  the 
e3q)erimental  systems. 

It  can  be  expected  that  as  engineers  become  a,yBxe  of  the 
numerous  apjd.icatlons  of  the  computer  in  their  respective  fields, 
they  will  begin  to  use  thmn  to  fhller  advantage.  Pszhaps  as  the 
novelty  of  the  equipment  wears  off  and  the  word  "computer"  ceases  to 
be  synonymous  with  myetexy  and  complexity  the  analog  computer  will 
coonand  widespread  use  in  engineering  organlntions  of  both  govemnent 
and  industry. 
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PRCBUauS  IN  TBE  AFTLICATICN  OF  lEAST  SQUARES 
TC  Fl’NCTICN  GUNERATOR  DESIGN 
By  John  G.  Tappert 

A  function  generator  is  e  circuit  or  network  used  in  an  analogue  cong;)utar 
to  store,  and  e^diibit  on  demand,  a  predetexrained  relationship  between  a  number 
of  variables.  Consider  the  relationship, 

f  (y,  s  ,  X  ....  X  )  =  0. 

1  £  n 

A  function  generator  used  to  solve  this  relationship  would  have  the  property 
of  exhibiting,  in  analogue  form,  the  value  of  y  which  satisfies  the  given 
equation  whenever  a  set  of  x's  eure  introduced  in  their  appropriate  analogue 
form.  Nomally,  it  is  required  that  the  relation^ip  be  satisfied  only  in  a 
limited  domain  and  within  a  specified  accuracy.  Further,  in  confuted  applica* 
tions,  the  functional  relationships  are  limited  to  those  which  are  single 
valued  and  have  finite  first  derivatives  of  the  dependent  variable. 

Function  generators  are  essential  to  many  military  conqputers.  Practical 
design  of  such  circuits  requires  that  efficient  use  be  made  of  available  com¬ 
puting  elesKints.  Cbe  of  the  most  sisqile,  yet  extremely  versatile  computing 
elements  is  the  fxmction  of  a  single  variable.  That  is,  an  arbitrary  function 
of  a  single  variable  is,  in  most  cases,  easily  and  reliably  instrumented  to  a 
high  degree  of  accuracy.  It  is  reasonable,  then,  in  designing  a  function 
generator  involving  many  variables  to  attempt  to  approximate  the  desired 
functional  relations  by  combinations  of  functions  of  single  variables.  Ihls 
paper  will  deal  with  the  mathematical  problems  Involved  in  the  determination 
of  optimum  functions  of  single  variables  in  higher  order  curve  fitting. 

The  techniques  for  applying  the  principle  of  least  squares  to  curve  fitting 
problems,  in  the  case  where  a  finite  niunber  of  parameters  are  to  be  pptimiaed, 
are  well  known.  If  the  residuals  con  be  expressed  as  linear  functions  of  the 
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parameters  to  be  deteimined,  then  these  parameters  can  be  determined  by  solring 
a  set  of  linear  equations*  In  our  present  problem*  the  unlmowns  are  functional 
relationships  rather  than  coefficients  of  known  ftmctions.  Howerer,  the  foxmer 
is  readily  reduced  to  the  latter  by  defining  the  functional  relationship  in 
terms  of  a  finite  number  of  parameters* 

Let  ns  consider  an  example  along  with  a  typical  set  of  conditions*  Let 
there  be  a  relation 

Yl  *  f  Xg.  x^).  (1) 

which  is  to  be  instrumented  in  the  form 

yg  =  (Xj^)  Fg  (Xg)  -J.  F^  (x^)*  (2) 

It  is  desired  to  determine  the  fimctions  F^*  Fg*  F^  idiich  minimize 

dXj  dXg  dx^  (3) 

oyer  some  specified  domain  in  X-sfiace*  However,  no  explicit  analytic  expression 
for  (1)  is  known]  but  it  is  possible*  at  considerable  cost  in  time  and  labor*  to 
determine  sets  of  data  which  satisfy  the  relation*  Further*  the  procedure 
enployed  in  deriving  these  sets  of  data  does  not  pexmit  direct  choice  of  apeelfie 
values  of  x's  or  y;  thus*  the  available  data  is  irregularly  spaced  in  all 
variables* 

It  might  seem  reasonable*  under  these  conditions*  that  if  we  took  a  repre¬ 
sentative  set  of  data  and  replaced  the  integral  (3)  by  a  summation*  the  problem 

could  be  solved*  But  suppose*  for  exanple*  that  among  the  selected  data  points 

no  two  values  of  say*  x^^  were  identical.  Then,  as  far  as  the  data  are  concerned, 
y^  could  unambiguously  be  represented  as  a  function  of  x^  alone*  and  a  solution 
with  all  the  residuals  vanishing  would  be  obtained  by  making  Fg  and  constants* 
Such  a  soluti(xi  obviously  would*  in  general*  not  satisfy  the  condition  (3)  and 
therefore  would  not  be  satisfactory* 
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These  oonditionet  then*  are  not  sufficient  to  solve  the  problem.  It  is 
necessary,  as  it  always  is  when  we  expect  to  define  continuous  functions  by  a 
finite  number  of  points,  to  assume  that  the  functions  have  that  rather  ill- 
defined  characteristic  of  being  ■well  behaved*.  Tor  our  present  puxpose.  we 
can  say  that  function  is  «well  behaved*  if  it  has  no  structural  peculiarities 
that  are  not  revealed  by  the  available  data.  This  property  of  *well-behavedness* 
is  akin  to  the  property  of  continuity  in  that  it  is  a  local  charaoteristic. 
involving  only  the  relation  of  data  points  in  a  finite  but  limited  region, 
ffe  are.  then,  faced  with  a  two-fold  problemt 

a.  To  define  a  funetlonal  relationship  in  terns  of  a  finite  number  of 
parameters. 

b*  To  define  a  functional  relationship  in  such  a  manner  that  it  is 

locally  *well-behaved*  at  evezy  point,  without  imposing  restrlotlons 
on  Its  gross  shape. 

There  are.  undoubtedly,  many  ways  in  which  these  two  conditions  can  be 
met  reasonably  well.  The  stethod  that  has  been  employed  successfully  at 
Ftankford  Arsenal  consists  of  defining  the  function  in  terms  of  a  finite 
number  of  cardinal  values  and  a  low  order  interpolation  law.  The  procedure 
in  defining  a  function  of  x  would  be  as  followsi 

Nonnalize  the  variable,  x.  so  that  its  excursion  is  bounded  by  convenient 
cardinal  numbers,  say.  0  to  10,  and  define  the  value  of  the  function  at  all 
the  cardinal  points  of  the  argument  by  a  parameter.  Then  the  function  would 
be  completely  defined  by 

f  (N  k)  -  2  On  'n* 

where  the  r's  are  the  cardinal  value  parameters,  the  a's  are  the  interpolation 
coefficients,  and  the  argument  is  written  as  a  cardinal  number.  N.  plus  a 
decimal  fraction,  k.  In  the  ease  of  Xagranglan  interpolation,  the  numerical 
values  for  the  a's  depend  on  k  only  and  their  assignment  as  coefficients  of 
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qpeolflc  F's  depends  only  on  N* 

4 

This  procedure  obriously  defines  a  functional  relation^ip  in  terms  of  a 
finite  number  of  parameters:  the  low  order  inteipolation  law  adequately  controls 
the  "well-beharednese*  of  the  function;  and  the  fact  that  with  low  order  inter¬ 
polation  all  3  coefficients  waniah  except  those  for  the  cardinal  values  in  tha 
imnedlate  neighborhood  of  the  argument  limits  the  restrictions  imposed  by  this 
definition  on  the  gross  shape  of  the  function* 

Ihe  problem  of  optimizing  functional  relationships  can*  thereforsi  be 
reduced  to  the  problem  of  optimizing  constants)  and  this  problem  in  turn  is 
amenable  to  solution  by  known  conventional  methods* 

X)q;>erlenee  has  shown  that  difficulties  can  arise  if  eertain  prseautions 
are  not  taken.  Some  of  these  are  discussed  below. 

It  is  very  desirable  to  define  the  mathematical  problem  in  sueh  e  manner 
that  the  solution  is  unique.  Consider  a  stnieturs 
y  • 

and  let  f..  and  t  be  functional  relations  which  result  in  a  sum  of  the  squares 
10  20 

of  the  residuals  which  no  other  funotioms  ean  reduce*  Then  there  are  a  whole 
set  of  functions t 


f 

2 


(X^)  s  k  f 

(V  •  i  f 


10 

20 


(X^). 

(X^). 


which  give  identically  the  same  residuals*  These  solutions  are*  of  course f 
identical  as  far  as  a  conq>uter  designsr  is  ccncexnsd*  but  are  different  solu¬ 
tions  mathematically*  and  therefore  leave  things  indeteailnate*  This  type  of 
difficulty  can  usually  be  eliminated  by  reducing  the  number  of  parameters  by 
one,  and  defining  the  function  at  one  of  the  cardinal  values  in  teau  of  anothert 


^1  ^^2^  *  ^1  ^^1^  ^  constant 
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This  fixes  the  scale  of  the  f^  function  without  restricting  its  shape,  and 
eliminates  the  multiplicity  of  equivalent  solutions.  In  more  complicated 
structures,  the  possibilities  of  such  multiple  solutions  may  be  less  obvious, 
and  their  elimination,  without  ioposing  unnecessaxy  restrictions,  can  call 
for  considerable  ingenuity. 

Every  effort  should  be  made  to  reduce  to  a  minimum  the  number  of  undeter¬ 
mined  parameters  required  to  define  a  function  generator  structure.  Ihls  not 
only  reduces  the  amovint  of  labor  required,  but  also  provides  the  best  insurance 
that  the  solution  is  not  unduly  influenced  by  the  particular  finite  selection 
of  data  which  are  taken  to  represent  a  bounded  continuum. 

A  designer,  in  setting  up  a  structure  to  be  investigated,  usually  has  a 
fair  idea  of  the  general  fom  of  the  functions  which  will  be  required  for  it. 
Say,  for  exanple,  that  there  is  reason  to  believe  that  one  of  the  functions 
will  turn  out  to  be  something  like  the  tangent  of  an  angle,  between  0^  and 
85^.  That  is, 

f  (z)  V  tan  z. 

To  adequately  represent  such  a  function  with  a  low  order  interpolation  law 
would  require  closely  spaced  pnzwmeter  values  in  the  region  of  large  values 
of  z,  but  could  tolerate  wide  spacing  in  the  region  of  small  values  of  x. 
Oonsiderable  economy  of  parameters  can  be  realized  by  a  change  of  variable, 

f  (x)  i  g  (tan  z). 

Since  the  g  function  would  be  alir.ost  linear,  it  presumably  could  be  adequately 
defined  by  a  small  nuid»er  of  parameters  and  still  retain  the  advantage  of  a 
low  order  intexpolation  law.  Alternative  ways  of  handling  this,  vdien  the 
change  in  variable  is  not  convenient,  amount  to  breaking  up  the  undetexmined 
function  into  two  terns  or  factors,  one  of  irtiieh  is  predetezmined.  Thus, 
there  is  no  loss  in  generality  in  defining 
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f  (z)  *  ten  X  •  gg  (x). 

Heret  again,  fewer  parametera  would  be  required  to  adequately  define  g^  or  £2 
than  would  be  required  to  define  f* 
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KAUBEtfiiriCIAirS  AS  m 
I|3r  Saul  Qom 


nnsDEucTion 

At  t«D  other  meetings  sponaored  1::|7  the  Office  of  Ordnance  Beeeareh, 
namely  the  ^posla  of  April  19^  and  iprll  199^  in  Hew  York  and  Chicago 
re^ctlvely,  there  were  talks  on  subjects  similar  to  the  one  I  am  called 
t:^n  to  discuss  today.  Since  Ik*.  lotkin,  idio  gave  the  address  in  Hew  Toric, 
ani  Professor  A.  A.  Bennett,  who  spoke  in  Chicago,  were  both,  like  myself, 
connected  with  the  Computing  laboratory  of  the  Ballistic  Besearch 
feboratorles  at  Aberdeen  Proving  Oround,  it  might  seem  that  they  have 
scd^nly  discovered  some  Mg  gap  in  their  diseuBBions  for  me  to  so^maid 
t^n. 


Ihis  they  de^jiitely  did  not  do. 

And  yet  1  am  sure  that  they  would  both  agree  that  the  automatic,  hi^ 
floeed,  digital  computers  have  had  such  a  dramatic  effect  on  both  oar 
ability  to  solve  our  problems  and  our  attitude  on  how  to  do  so  that  this 
effect  will  bear  describing  again  and  again. 

However,  you  will  understand  why,  in  this  description,  I  cAiooas 
different  examples  from  theirs. 

Both  referred  to  the  variety  of  machines,  including  Aiiac»  Kvbc, 
Orivac,  etc.,  and  Dr.  lotkin  desc^bed  their  main  characteristics. 

Professor  Bennett  described  the  variety  of  application  and  Dr.  lotkin 
gave  their  mathematical  categories.  Tinally,  PTofessor  Bennett  gave  the 
percentage  of  computing  time  dented  to  ead:  class  of  application  and 
spoke  about  the  effect  of  hi^  speed  computers  on  matbsmaties  itself. 

I  want  to  touch  tqion  the  sources  and  nature  of  the  requests  for  oar 
services. 

Within  B(BL  we  get  problems  from  the  Drterior  Ballistics  laboratory. 
The  Interior  Ballistics  Zabozatory,  The  Ballistic  Nsasoramants  Xaboxatoxy, 
The  Terminal  Ballistics  laboratory.  The  Weapon  Sbrstems  Xaboxmtoxy,  and 
problems  from  the  Computing  laboratory  itself. 

Almost  one  third  of  the  projects  are  concerned  with  eospotations  of 
firing  and  bombing  tables.  Frm  exterior  ballistics  come  many  data 
reduction  problens  Involvltg  smoothing,  least  square  fitting,  and  the 
like,  as  well  as  a  number  of  air  flow  problems.  Interior  ballistics 
has  been  yielding  problems  in  the  determination  of  trajectories,  burning 
rates,  and  molecular  structures.  Beside  snaerous  data  redaction  proMsms, 
ballistic  measurements  give  rise  to  various  theoretical  problsms  conoexnsd 
with  measuring  systems.  Ws  will  see  an  example  in  a  aioment.  Berminal 
ballistic  problems  give  rise  to  a  host  of  theoretiml  stadias  involvii^ 
elasticity,  peoaetrabillty,  fbogoentatlon,  and  distribution  of  fkagnents. 
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Ihe  pro'blems  of  weapon  systems  erml’jation  are  largely. concemed  with 
pro'ba'bility,  statistics,  and  game  theory.  linally,  in  addition  to  the 
firing  and  homhing  tables,  The  Cospating  laboratory  has  many  problems, 
mathematiced  and  logical,  not  directly  tied  to  specific  computations, 
concerned  vdth  efficient  use  of  the  hi^^  speed  computers.  More  will  be 
said  on  this  subject  later. 

The  sources  of  problems  outside  SBL  are  the  arsenals,  other  Gbsement 
Installations,  Qovemment  contractors,  and  universities,  ^om  the  arsenala 
we  have  been  getting  large  data  reduction  problems,  often  by  comlicated 
least  squares  procedures,  and  computations  of  complicated,  trajectories; 
other  problems  fron  Qovemment  installations  and  Qovemment  contractors 
have  involved  computations  of  comlicated  multiple  integTsds,  hi^ 
precision  solutions  of  complicated  numerical  equations,  bo^mdary  value 
problems,  and  solutions  of  median  sized  systems  of  linear  equations. 

Problems  fiom  'Ture  mathematical  research  come  to  us  by  way  of  research 
contracts  awarded  by  the  Office  of  Ordnance  Research  itself.  An 
example  of  this  type  will  also  be  touched  upon  shortly. 
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An  example  of  how  an  advanced  mathematical  viewpoint  can  aid  in  the 
choice  of  an  economical  computing  method  by  examining  the  variables  in 
a  range  in  iMch  they  are  x^sic^Iy  meaningless  is  the  following.  The 
interior  balllsticians,  in  studying  the  motion  of  a  projectile  within  a 
gm  after  the  dmurge  has  bozned  out ,  have  a  phgreiflal  aodal  laadinc  to 
^  tzajaotoiy  aquation: 


where  Y  la  a  constant  in  the  nei^bozhood  of  I.3  and  Xand  U  are 
dimensionless  quantities  related  respectively  to  the  free  volone  behind 
the  projectile  and  the  pzpjectile  velocity.  It  is  desirable  to  have  a 
table  of  solution  curves  for  various  seta  of  initial  conditions 
corresponding  to  different  guns,  different  chargee,  different  projectile 
weii^ts,  etc.  Aram  the  vezy  statement  of  the  problem  it  is  evident  that 
the  equation  mnst  be  physl^Iy  meaninglees  in  a  xwighbozhood  of  x.*  O  • 
and,  indeed,  the  computed  values  are  only  of  physical  interest  in  the 


zaage 
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Bevartheless,  to  chooae  an  efficient  organisation  of  Ihe  computation  calls 
for  a  view  "in  the  large"  of  all  the  solutions  in  the  x,  a  -  plane,  whose 
diatribution  la  greatly  influenced  by  the  nature  of  the  singularity  at 
the  origin.  Sy  an  ovezall  analysis  of  the  type  much  used  in  nonr.linear 
medianicf  one  can  quickly  akmt^  the  feunlly  of  solutions  from  the  nature 
of  the  singularity,  the  aero  isocline,  the  infinite  iaocllne,  and  the 
locos  of  inflections;  all  these  may  ba  dstszminsd  directly  fxm  the 
emtion.  Ihs  situation  is  roughly  acasnarisod  in  the  accoatpanylng 
diagBan: 


ktleiaat 
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Because  of  this  view  of  the  solutions,  it  was  decided  to  organize  the 
nuraerical  solution  by  machine  by  1.  •  beginning  each  solution  at 

2  V  • 

u  ■  X  '  ,  proceeding  to  the  right  until  u  ■  >1.  2.  •  Sultching  variables 

and  continuing  on  to  the  edge  of  the  deslrxl  region,  3.  -  proceeding  left 

from  the  same  initial  point  in  the  same  way  swi telling  when  u  ■  1*  This 
procedure  was  carried  out  automatically. 

A  BAIUSTIC  MEASUREMENTS 

The  physicists  concerned  with  balllstlo  measurements  somatines  measure 
a  distance  bj'  sending  out  a  radio  signal,  frequency-modulated  hy  a  number 
of  different  tones.  Ihe  signal  is  automatically  returned  from  the  object 
whose  dintance  is  being  measured,  and  the  phase  differences  in  the 
different  modulating  frequencies,  also  called  nui^  channels,  are  used 
like  ttie  dial  readings  on  a  gas  meter  to  give,  so  to  speak,  successive 
decimal  digits  in  the  distance.  When,  however,  there  are  si^perfluous 
return  signals,  due  to  such  things  as  uiwanted  reflections,  a  distortion 
of  these  range  channels  is  the  result,  causing  an  error  in  the  estimate 
of  distance.  To  get  an  idea  of  the  magnitude  of  such  errors,  a  simi^fled 
model  with  one  unwanted  reflection  was  assumed,  and  the  reflected  signal 
was  supposed  to  have  a  strength  as  laige  as  that  of  the  direct  signal 
itself. 


The  Coasting  Laboratozy  had,  then,  to  produce  accurate  answers  to 
the  foUowlngi 

Find  the  phase  shifts,  ^  aad^,  for  two  values  of  the 
deviation  ratio,  p,  and  elf^t  values  saoh  of  9^,  v,  and 
^  ,  having  the  following  meanings  - 


-  Oq  ■  carrier  phase  lag  at  the  receiver  of  the  reflestsd 

wave  relative  to  the  direct  wave, 

V  -  a  quantity  proportional  to  the  path  difference 

^  -  a  quantity  connected  with  phase  lag  between  modulating 
^  frequencies. 


where  ^  and  h(2  are  given  hjr: 
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<^,=  -^S*rt  a  V  (H., -M.)  ^  (S, i. v) , 

H,  *  -~  ^-■-  i  ( X  -X)  OS  (3i  i-  v)  ^ 

Q^=^s^aov  (X-rX)  (/o3,-  -f-zov-J^ 


=  tz^lS^y  +/-2.  (x  ~X>  Crs  (/i'Z,-  V' vl 
y 


i*/ 


'oijXi )  6  Sre  SueX  as  "/i^ieAJ 

»  (iJr,z-/)7r 
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(Diere  are  essentially  nunerical  equations,  althou^  the  inter¬ 
mediate  variables  must  be  integers.  For  different  assigned 
integers  for  M ,  the  last  equation  most  be  solved  for  B  to  ^eld 
a  corresponding  Ihe  main  pro'U.em  is  essentially  one  of  or^niz- 
ing  the  eunangement  of  solution  of  somevdiere  in  the  order  of  10,000 
such  equations  in  such  a  vay  that  the  answers  are  properly  grou^ 
to  find  0]^,  Hi,  Og,  and  H2  and  hence 0(1  and<^2*  combimtlon  of 

bookeeping  ana  txaffic  control  with  (xmputed  intemediate  results  is 
one  of  the  main  capabilities  of  the  automatic  speed  machines. 

Ihe  numerical  methMs  adopted  would  have  taken  a  skilled  and  tireless 
hand  computer  at  least  one  man  year.  Less  than  thirty  hours  of 
time  was  actually  Involved. 

A  STATISTICAL  FBCBiat 

In  evaluating  the  effectiveness  of  a  gun  and  ammunition  it  is 
customaxy  to  use  the  standard  deviation  ftxn  the  mawn  point  of  impact, 
both  in  tange  and  deflection.  Saaqole  variances  or  standard  deviations 
are  customarily  used  to  estimate  the  corresponding  population  pEurasteters. 
However,  in  order  to  reduce  the  effect  during  xatge  firing  of  shifts  in 
the  mean  point  of  impact  due  to  wind,  A.  A.  Bennett  around  1919  pxvposed 
that  successive  differences  in  the  observed  ranges  and  deflections  be  used 
to  obtain  estimates  of  the  aforementioned  parameters.  This  haa  led  to  the 
problem  of  computing  the  probability  distribution  of  such  suople  ataiiatios 
aa: 


,  where 


,  aiid 


lifecif,  ,pf  cr-urao. 


A^X.  =  A'-'X 

and  th'i  are  Ir.  .ev!i,'.dent  random  sa,vples  fro.m  a  normal  populrtiou. 


It  ie  desirable  to  have  tables  of  these  distributions  for  somple  sizes, 
n,  from  3  to  29,  and  order  of  difference,  p,  from  1  to  9.  ^  distributions 
ttoselves  vmre  completely  characterized  by  von  Heomann  as  long  ago  ae  19^1; 
the  problem  ie  to  dioose  a  reasonable  methM  of  computing  them. 
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At  least  nine  methods,  ranging  f^om  the  nalve-hut-not-to-t>e-i^ored 
to  the  sophisticated,  haw  been  considered.  The  method  being  inyestigated 
at  present  is  somewhat  typical  in  that  it  probably  would  not  have  been 
considered  if  the  hi^-speed,  aatomatie  machines  had  not  been  awilable. 

First  an  orthogonal  transformation  is  applied  to  reduce  the  probleB 
to  that  of  computing  the  distribution  of 


?// 


The  diaraeteristic  values,  Bj,,  for  the  various  n  and  p  have  bean  computed. 

The  method  under  ooxxsidexatlon  begins  with  von  Ssunann' s  relation  for 
the  probability  density,  T  0) ,  applicable  idiea  m  is  even: 


Bus  a  dBstoh  of  f  (p)  end  ita  l/an-l  derimtlves,  ^tm  a *6,  wmdd  look 
as  followa: 


'Corn 
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If  one  expands  the  non-aero  expression  for  the  -y—  -  I'st  derivative 
of  F  (p)  into  fractional  power  series  about  each  of  tne  singular  points, 

the  radius  of  convergence  will  be  the  smaller  of  -  B^_j  and  Vl  ■  ®v. 

Integrating  these  series  tern  hy  term  should  permit  us  to  advance  along  the 
p-axis  by  a  method  of  analytic  extension  applied  to  all  the  derivatives 
simultaneously  down  to  P.  In  all,  this  will  call  for  the  use  of  almost 
7^00  power  series  and  almost  as  many  polynomials,  which  should  be  handled 
automatically  in  interrelated  groups  to  economize  on  the  requiremants  for 
simultaneous  storage  of  coefficients. 

A  PROBUSM  FBOH  MATHEMATICAL  tfiSEARCH 


The  following  problem  ooneemed  with  the  theory  of  algebraic  numbers 
came  to  us  from  the  Office  of  Ordnance  Research  in  support  of  a  oontraet 
it  had  awarded.  A  nuiriber^  is  called  a  cubio  algebraic  number  if  it 
satisfies  an  equation  of  we  forms 


solutions. 


the 


conjugates  of 


called  the  trace  of 
^je^^^called  the  norm  of  ^ 
t,  T,  and  N  rational  niudsers,  and 


satisfies  no  such  equation  of  lower  dsgroe. 


If  t,  \  and  N  are  also  Integers,  Is  an  Integral  algebraic  number 
If,  in  addition  N  ■  jt  1,  ^  Is  called  a  unit.  Like  the  mudbers  ^  1, 
units  divide  exactly  the  other  algebraic  Integers  of  their  domain* 


Professor  H*  Cohn  of  Wayne  Ikiiverslty  designed  an  algorithm 
which  within  a  cubic  field  produces  mangr  algebraio  Integers  with  small 
norms,  and  must  yield  at  least  one  unit.  It  begins  with  a  suitablj'  chosen 
three  bgr  four  matrix  of  algebraic  nunibere 
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with  sign  pattern 


and  P  +  Q  +  R  +  S»0 

One  Cb.i  choose  a  column,  add  It  to  a  second,  and  subtract  It  fron  a  third* 
Of  the  twenty  four  possible  results,  between  three  and  six  will  Maintain 
the  sign  pattern,  and  the  algorithm  consists  In  branching  out  from  these 
by  the  same  process.  However  a  matrix  Is  the  end  of  Its  line  if  some 
suitable  rearrangement  of  its  rows  and  columns  makes  Its  rows  individually 
proportional  to  those  of  a  predecessor.  Ihe  faoton  of  proportionality 
are  units* 

In  order  to  make  more  modest  demands  on  storage  the  sponaor  replaced 
this  branching  algorithm  by  one  whose  results  appeared  in  a  linear  order* 
The  new  algorithm  was  to  be  applied  to  the  so-called  *(qrolio  eiibie 
fields  with  prime  conductor*,  beginning  with  aatrlees  siivly  related  to 
the  solutions  oft 


a  prime  intofuv 


irtwre 


/  and 

V  ylX  , 

The  unique  solutions  of  /7  y*  ^  ^ 

It  was  found  possible,  using  about  e  thouaand  atorage  poeitlona  on 
one  of  our  maoUnes  to  cany  out  all  operatlona  oxaetly  oxoopt  for  thoM 
required  to  verify  the  sign  patterns  of  the  matrloes*  If  the  maohiaee 
had  more  Internal  storage  available,  ^ee  sign  dlsorlxliaitlona  ee«d4 
aleo  have  bem  made  ind^ndeot  of  approxlmaUcne*  Ae  it  was,  automatie 
error  controls  were  carried  along  with  the  computations* 


The  fact  that  the  essential  oosputatlons  coiQd  be  oarried  out  to 
what  is  effectively  infinite  precision  to  produoe  some  3^00  algebrajn 
Integers  together  with  their  determining  equationB  for  all  prliias  J[ 
up  to  $00  was  made  possible  by  the  use  of  ^  so-ealled  *stiuet«rs 
constants '  t 


f  f  sttoh  that 

t/,U,  tA  , 
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nftfiirO  y : 


By  use  of  thef  e,  every  int,9ger  forr«d  and  its  defining  equation  could 
be  given  as  follov/sj 

c)—  ^3 

a,  b,  and  c  lational  integers, 

t  -  (a  , 


where 


,<r,= 


This  problem  thus  had  the  Interesting  and  paradoscical  feature  of 
including  both  infinite  precision  conpitations  and  automatie  ezror 
analyses , 

THE  GENEHAL  PROBLEM  OF  ORGANIZING  MACHINE  USE  EFFICIENTLI 


The  problems  described  so  far,  although  diverse  in  eoerce  and  beck* 
ground,  have  in  common  a  kind  of  dynamic  logic  ve  have  already  referred 
to  as  boc^keeping  and  traffic  control.  They  have  all  included  the 
mechanization  of  decisions  on  ho\i  to  continue,  where  the  criteria  are,  of 
coxirse,  fixed.  In  this  sense  the  machines  can  think,  and  the  thinking 
can  be  on  a  large  scale  and  very  intricate,  but  must  be  routine  as 
opposed  to  creative  or  executive.  And  in  most  problems,  even  though 
tliey  are  computational,  a  large  part  of  the  work  is  this  making  of  routine 
decisions.  This  is  even  true  when  the  problem  is  not,  strictly  speaking, 
computational,  as,  for  Instance,  when  we  set  the  machines  up  to  jday 
war  games.  It  is  perhaps  for  this  reason  that  programs  for  the  machines 
are  commonly  called  routines.  One  of  the  most  surprising  facts  to 

emerge  from  constant  attempts  to  extend  the  range  of  applicability  of 
these  machines  is  that  vast  portions  of  our  judgments  are  routine  and 
can  therefore  be  mechanized.  Perhaps  these  machines  should  not  be 
called  conpiuters  at  all|  they  can  do  general  information  transforming 

of  which  there  are  other  forms  than  the  numerical  one  of  computing} 
translating  languages  is  an  example. 

An  example  of  a  large  scale  information  transforming  job  already 
mentioned  as  a  primary  function  of  BRL’s  Computing  Laboratory  is  the 
following:  Given  large  quantities  of  range  data,  produce  a  firing 
table,  whol?  procedure  can  be  thought  of  as  one  huge  routine  wb.ose 

subreutinos  at  present  use  different  machines  and  groups  of  people. 

The  poqile  are  used  less  and  less  within  this  routine  and  more  and  more 
in  plaij;lng  its  mechanization.  For  example,  a  big  part  of  the  hand 
typing  of  the  tables  haa  been  eliminated  because  the  extraction  of 
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the  tabulated  ififoxtsation  le  perfoaned  aectaaaioally  upon  the  oonpated 
tmjeetoriee  and  prepared  meehanloally  fbr  eatoaatle  tgrpii^S* 

■aoUnee  with  larfie  internal  etom^e  become  etfailable  to  oe,  eery 
little  of  the  table  production  pxooednxe  will  be  left  umieehaals^; 
eren  tihecking  for  bad  data  can  be  eede  routine. 

At  present,  then,  more  time  is  epent  preparing  probleme  for 
computation  than  it  tataes  to  do  the  eonpatatioa.  9ie  problene  are 
literally  quidoer  done  -than  said.  Nonthe  are  needed  to  prepare, 
against  hours  to  compute.  Qiis  ie  becauee  aneh  judgeent  le  required 
in  the  preparation  of  routinee. 

And  yet,  as  usual,  a  large  portion  of  this  Judgeent  is  itself 
routine,  end  can  again  be  aedhaniaad.  Ws  diould  be  able  to  mliat  the 
■adiines  to  help  in  preparing  their  own  probleme  for  their  own  eolntion, 
end  thne  poesiuy  reduoe  this  preptffatien  time  from  months  to  mlnntes. 

Let  ue  pease,  then,  to  break  down  this  preparation  tiae.  bhea  a 
problem  ie  reoeirod,  it  must  go  throu^d^  ^enr  staips  before  all  tha 
actual  oooqptttations  are  perfMMd.  fteea  may  be  eallod  analyela, 
error  analyela,  oodlng,  and  eede  cheddag. 

Cie  eaalyeie  etegs  inoludee  chooelng  the  methods  to  bo  need,  l.o. 
ohoosing  the  aeda  subrMtlnee  la  the  routine  together  with  their  logieel 
ordering,  and  doing  such  analytic  things  ae  say,  substituting  one  power 
eeries  into  another.  An  example  of  the  first  of  analysis  is  the 
posing  of  a  dlfferenoe  aquation  approxinatioa  to  a  differential  equation 
(the  only  etrietly  nomericel  opemtlone  theeo  ■Mhimee  aerdsni  are  adAltiem, 
mbteaetlom,  aultSiaiaation,  and  dirialea) .  Mm  immIV  aelle  M  sisMtlmm 
Judgsment  end  would  rarely  be  meofaealesd. 

Bm  second  stage,  error  aaalysie,  al«>  has  a  large  element  of  remtine 
thinking  in  It.  It  dionlA  therefore  be  possible  in  many  oases  to  erremge 
for  the  error  anelysie  end  oontrol  to  be  ourried  on  al^  with  the 
oomputatloa  by  the  merhine.  \fe  hare  mentioned  this  eort  of  thii^  in 
ooaaeotion  with  the  problem  fkom  mathsmatlaal  research. 

lha  coding  etage,  ae  everyone  is  aware,  alee  contains  a  lAxge  eUment 
of  the  routine  elerioal  type  of  thood^t,  eapoelelly  if  well  known  aetbida 
are  being  applied.  M  routine  element  here  eoaoiste  in  putting  together 
the  iaetrocttoae  themaelves.  But  our  mechlnes  etore  their  orders 
exactly  as  they  store  their  nasbero,  so  that  routina  operations  on  their 
own  orders  can  be  carried  out  modiBnically.  fhio  na\vi  it  poeeible, 
e^iin  if  the  storage  it  suffioiently  extended,  to  hate  ttie  nadilaee 
•tttomatioally  aeomble  their  om  routines  from  previoualy  coded  sni^ 
routinee.  Oombining  this  poesibility  with  the  aadiinee*  oepability  of 
doing  routine  translation,  this  means  that  It  is  possible  in  a  multi- 
machine  instaUation  each  as  SL*  s  Ooqpating  Xnbormtory  to  code  ptoblemo 
in  a  nnivoreal  abbreviated  code,  independent  of  the  purtioulnr  nmohiae. 
fte  nadilnes  thamoelves  would  aeemBble  their  oedes  itm  previouely  eeied 


Gonferenca  of  Arsenal  Mathematicians 


51 


sahroatises,  and  then  translate  the  routines  into  their  own  code,  naking 
their  own  routine  choices  of  storage  assignments.  This  is  a  unimal, 
seni-euitomatic  coding,  kxperlnents  are  being  run  on  the  spare  time  left 
1y  production  commitments  to  see  how  expensive  sudi  a  procedure  would  he 
in  both  storage  and  time. 

Ihe  last  stage,  code-checking,  can  take  considerable  time  sinoe 
hcoan  errors  in  coding  are  almost  teund  to  appear  in  long  problems, 
sudi  as  some  we  have  described.  The  existence  of  a  large  library  of 
standard  subroutines  would  do  much  to  alleviate  this  situation;  the 
existence  of  universal,  semi-automatic  coding  procedure  would  do 
even  more . 

With  these  considerations  in  mind,  it  is,  I  believe,  apparent  why  so 
much  has  been  aaid  about  the  effect  these  madilnes  may  have  on  the  d^lq^ 
ment  of  mathanatles  and  on  the  type  of  personnel  needed  to  pxogrsm  them. 
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COMPUTING  METHODS  TOR  THE  SOITTICN  CT 
SrSTElG  OF  DIFFERSRTIAL  EQUATIONS 
F.  J.  Murray 

The  Office  of  Ordnance  Beaearch.  has ‘kindly  giran  me  an  opportzinity  to 
discuss  the  work  which  I  have  done  in  connection  with  conoputing  procedures 
for  differential  equationst  and  I  am  grateful  for  this  invitation* 

The  experience  of  any  individual  is  necessarily  limited,  and  the  methods 
which  have  interested  me  have  been  developed  in  connection  with  a  certain 
range  of  problems  and  certain  types  of  confuting  facilities.  I  must  apologize 
then  if  this  range  is  not  coincident  with  your  range  of  interest. 

Most  of  this  work  arose  as  a  consequence  of  consultation  with  the  United 
States  Navy  Special  Devices  Center,  and  certain  theoretical  developments  were 
carried  on  later  under  the  auspices  of  the  Office  of  Naval  Research.  Special 
Devices  has  had  a  remarkable  and  interestixzg  range  of  problems*  It  has  had 
considerable  responsibility  for  the  development  of  design  conputational  procedorea 
for  air  vehicles  and  airborne  systems.  As  you  are  aware,  the  original  RBAC 
equipment  was  developed  under  SDC  auspices. 

The  application  of  this  equipment  is  dependent  upon  the  analysis  of  many 
mathematical  problems,  such  as  the  requirements  of  real  time  simulation,  overall 
system  analysis,  and  the  use  of  digital  cheek  solutions.  The  Cyclone  Labora¬ 
tories  of  the  Reeves  Instrument  Coxporation  pioneered  the  use  of  automatic 
digital  eooqputing  equipment  for  check  solutions  for  analog  ccsqaitatiaii.  The 
experience  with  the  digital  check  solutions  for  analog  equipment  was  later 
useful  for  developments  such  as  the  present  project  for  the  development  of  an- 
operational  flight  trainer  using  digital  coqputation* 

In  certain  of  these  cases  the  actual  problem  in  applied  mathematics  is 
more  general  than  the  specific  solution  of  an  individual  system  of  differential 
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equatlona*  The  problems  involve  classes  of  systems*  l^>r  instance •  the  develop¬ 
ment  of  analog  con5)uting  equiprent  involves  questions  concerning  ranges  of 
behavior,  and,  similarly,  the  major  objective  in  the  development  of  an  oper¬ 
ational  flflght  trainer  using  digital  computation  is  a  flexibility  idxich  will 
permit  the  consideration  of  a  reasonable  range  of  aircraft*  Ihus,  it  is  char¬ 
acteristic  of  these  problems  that  while  they  arc  not  spoeifically  tied  to 
individual  systems,  nevertheless,  like  almost  all  problems  in  applied  mathe¬ 
matics,  they  do  have  precise  limits* 

It  seems  to  me  that  the  first  responsibility  of  an  applied  mathematician 
is  to  solve  the  given  problem  with  precisely  the  generality  given  to  him* 

If  more  than  one  solution  is  possible,  his  objective  is  to  obtain  the  optimum 
solution.  Until  the  given  problem  is  solved,  his  major  interest  should  not 
become  involved  in  any  generalixation  of  the  given  problem*  After  a  specific 
problem  of  a  certain  type  has  been  solved,  the  solutioa  itself  may  give  same 
light  on  effective  generalizations  and  tie-ins*  The  solution  may  be  a  small 
window  which  will  permit  us  to  see  generalizations  that  are  really  valuable* 

Ihe  work  for  SDC  indicated  above  has  led  to  two  areas  of  developiMnt 
which  I  e^qplored  under  CNB  auspices  with  two  colleagues*  Dr.  Killer  and  I 
have  developed  a  general  error  theory  for  differential  analyzers*  (Cf*  Ref* 

)  This  is  adequate  for  digital  procedures  as  well  as  the  more 
difficult  continuous  coirputer*  J*  Winson  has  modified  this  theory  euid  applied 
the  result  to  "linear*  MIC  equipment  with  very  interesting  results*  (Ref*  ) 

4 

Other  developments  of  this  theory  were  undertaken  by  the  Advisory  Board  on 
Siaulation  of  the  Universl'ty  of  Chiceigo*  The  latter  will  also  publish  sn 
extended  version  of  the  theory*  Dr*  Brock  and  I  have  worked  on  certain 
numerical  eooputing  procedures  idiich  are  applicable  to  the  types  of  systems  of 
differential  equations  fdiich  eppear  in  problems  for  air  vehicles*  (Hsf*|g  and§J) 
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It  present,  it  aeems  to  me  that  it  is  less  expensire  and  more  eanvenient 
to  study  a  system  of  differential  equations  of  the  type  ehich  occurs  in  air 
rehicle  problens  hy  the  use  of  a  continuous  differential  analyzer  properly 
aippleeieBted  by  digital  computation,  than  any  olher  method*  Parameters  may 
be  varied,  equations  'simplified*,  suid  the  equlyalent  of  major  design  changes 
carried  out  with  ease,  The  person  with  engineering  responsibility  for  the 
design  is  at  all  times  in  intimate  contact  with  his  problem*  TUrthexmore, 
real  time  simulation  at  present  seems  to  require  the  use  of  continuous  equip¬ 
ment*  I  was  pleased  to  note  Ifr.  Wood  in  the  present  meeting  has  indq>endently 
made  aimllar  remarks  concenilng  the  applicability  of  continuous  coiqputers* 

There  eoe,  however,  many  problems  associated  with  the  use  of  such  equipment* 
To  check  the  initial  setup  and  at  the  same  time  obtain  a  prelimlneury  test  of 
accuracy,  a  digital  check  problem  is  very  useful*  The  size  of  these  problems 
maksa  aatomatlc  computen  neee— iry,  and  eritiaal  reamitm  ana  be  ehaahad  bgr  a 
digital  ccnputation* 

Ihlle  the  use  of  digital  cheek  computation  is  probably  the  best  manner  of 
detemlnlng  the  accuracy  of  the  results  of  a  specific  study,  it  is  also  desirable 
to  obtain  a  general  theory  for  the  accuracy  of  continuous  differential  analyzers* 
Mathematically,  this  means  that  we  must  obtain  a  conparison  between  the  solu¬ 
tions  of  two  sjratems  of  differential  equations  which  in  soaie  sense  approximate 
each  other* 

Snob  a  comparison  is  desirable  and  essential  in  many  other  problems  besides 
error  studies  for  differential  analyzers*  .In  almost  every  physical  applica¬ 
tion  of  ordinary  differential  equations,  the  corresponding  prablem  arises  of 
comparing  the  solutions  of  a  known  system  of  differential  equations  with  the 
solutions  of  another  system  of  differential  equations  vdiich  approximates  the 
first  system* 
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Now  if  this  problem  involved  only  differential  equations  of  the  same  order* 
the  usual  existence  theory  of  ordinary  differential  equations  which  discusses 
the  dependence  on  parameters  of  two  systems  of  the  same  order  would  give  an 
adequate  answer*  But  both  in  the  ease  of  the  error  analysis  and  the  other 
applications  referred  to,  one  has  the  situation  that  the  two  systems  to  be 
compared  are  of  different  order.  Thus,  mathematically,  one  has  the  problem 
of  comparing  two  systems  of  differential  equations 

-  0  Oj  ■  0 

where  the  second  system  may  be  of  higher  order  than  the  first,  but  the  extra 
hl^er  derivatives  appear  in  the  second  system  with  relatively  small  partial 
derivatives.  This  is  the  problem  which  Dr.  K.  S.  Miller  and  I  solved  in  our 
joint  paper  in  the  M.I.T*  Journal  of  Mathematics  and  Physics.  (Ref.  ) 

The  review  in  Mathwwwtlcal  Reviews  of  this  paper  gives  the  incorrect  ijopression 
that  only  a  much  more  specialized  problem  was  considered.  Since  the  major 
objective  of  the  paper  was  exactly  thie  mathematical  problen,  it  was  quite 
disheartening  to  find  that  the  generality  of  the  result  was  ignored* 

In  order  to  tackle  this  problem  associated  with  an  increase  in  order, 
one  must  introduce  a  new  technique.  The  basis  of  this  technique  is  the 
followingi 

Around  a  given  fixed  solution,  one  may  "linearize*  the  given  system  of 
equations,  using  the  classical  existence  theory  procedure.  "Linearization* 
of  course  refers  to  a  process  of  replacing  the  given  system  by  a  linear  aystmn 
in  idilch  one  has  constant  coefficients*  Now  the  usual  rough  method  of  attacldag 
this  is  the  followings  Che  starts  with  a  single  solution  given  over  a  fixed 
interval  of  the  independent  variable  t*  Che  is  trying  to  study  the  solutions 
with  initial  conditions  close  to  this  one*  The  fixed  t  interval  is  divided 
into  shorter  intervals  and  on  each  of  these  the  given  system  is  replseed  by 
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a  aet  of  linear  equatlona  with  constant  coeffieients  and  appropriate  forcing 
tema*  It  ia  not  vmeooimon  to  use  this  system  to  study  the  original*  the  claim 
then  being  that  if  one  repeatedly  subdivided  the  given  interval  and  took  Initial 
conditions  very  near  to  those  of  the  given  solution*  one  would  (q;>prozimate  the 
various  other  solutions. 

It  does  not  seem  to  be  widely  appreciated  that  the  classical  theory  of 
differential  equations  permits  one  to  "linearize*  in  a  much  more  effective 
manner.  He  may  have  to  use  a  subdivision  of  the  t  interval  into  a  number  of 
saaaller  Intervals  and  use  as  the  basic  approximation  a  niunber  of  linear  eystems 
with  constant  coefficients,  but  one  can  base  the  discussion  on  a  fixed  aub> 
division  of  the  t  interval  and  the  range  of  initial  conditions  can  also  be 
considered  to  be  fixed  and  of  finite  extent. 

Briefly*  the  process  consists  of  introducing  a  system  of  equations* 

*  0* 

with  a  parameter  *>wltb  the  following  propertiest  Ihe  aystma* 

F^(..*  *0)  •  0* 

is  the  besie  approximation*  i.e.  it  is  a  linear  ss^stem  with  coefficients  idiich 
are  constant  on  intervals  and  change  at  most  a  finite  nuiid>er  of  times  on  the 
given  interval.  Ihere  is  also  a  value  (k^for  the  parameter  such  that 

FjC...  0 

is  the  original*  possibly  non-linear,  eystem  that  we  wiah  to  study, 
fe  now  consider  the  system 

Fj(**.  *s(9  *  0 

Cxls  is  a  system  of  differential  equations  on  unknown  functions  x^(t*0L)  which 
normally  can  be  considered  to  be  analytic  in  the  parameter  CL  A  solution 

fXinctlon  x^(t,(C.)  then  can  be  expanded  ia  jiowers  of  CL  arovind  (t*  0.  Thus 

Xi(t*a)  «  Xi(t*  0)  (t*  o)a  (t*  o)ct^  ^ 


•  •  •  t 
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Tbe  setup  mist  be  such  that  this  aeries  ecuTerges  for  dP  This  is  the 
single  assuoption  made  in  the  lAole  process  and  this  question  has  to  be  explored 
in  each  indiridual  case*  HoseTer*  the  theory  shows  that  eren  if  one  used  a 
linearization  with  coefficients  constant  for  the  ehole  interralt  there  is  a 
range  of  initial  conditions  for  which  one  has  conwergence.  It  nay*  howerer* 
be  more  convenient  and  siffiificant  to  linearize  on  subinterrals* 

The  above  discussion  is  based  on  the  existence  of  a  single  ■clothes  line* 
solution  along  which  the  regions  are  pinned.  X.  S.  Killer  and  I  have  written 
a  book  on  Xzistenee  Theorems  idilch  will  appear  very  shortly.  (Cf.  Ref* 

Ih  this  we  indicate  that  there  are  eesentially  two  classical  situations  for 
exlstenee  theoresis*  one  of  these  is  the  "initial  neighborhood*  or  *ln  the 
■sail*  theory  t  the  other  requires  the  existence  of  such  a  running  solution* 

In  practice  the  existence  of  the  running  solution  or  a  number  of  such  can 
be  established  by  caaputation«  and  thus  cosputation  and  theory  ean  be  utilised 
to  give  a  cooplete  theory  which  is  inaccessible  by  theory  alone*  ConTeraely* 
the  theory  givee  preolee  local  information  not  available  from  cooputatlon* 

This  book  is  being  publiabed  by  tbe  Hew  York  University  Press* 


How  the  linearized  equation  gives  tbe  critical  information  required  for 


the  linearization*  The  x^(ttdL)  satisfies  the  linearized  system  while  the 
satisfies  a  linear  system  of  equations  with  the  same  homogeneou  part 
but  different  forcing  terms*  and  so  do  the  and  the  bicker  partials. 


ThnSf  our  linearization  process  refers  the  original  problem  to  the  study 
of  tbe  x^t  *  lach  of  &ese  sets  of  functions  satisfies  a  linear 
egretem  of  equations  with  constant  or  interval-constant  coefficients*  The 


systems  for  the  different  sets  differ  only  in  the  forcing  terms* 


How  the  above  linearization  permits  us  to  attack  ortr  original  problem 
of  cosyaring  the  solutions  of  different  sj^stams  by  linearizing  both  the  glvwi 
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ejrbtem  and  the  approximating  eystem*  Aa  far  aa  basic  stability  is  eoneemed«  'fte 
Tariatlon  in  the  forcing  terns  is  insignificant*  Thus*  the  final  result  of 
the  linearization  process  is  to  set  up'  a  corresponding  problem  on  a  linear 
system  of  differential  eq.uatlon8  with  constant  coefficients  or  interral-eonstaat 
coefficients* 

This  last  problem  can  be  directly  attacked*  In  a  linear  aystemt  the 
increase  in  order  with  snail  coefficients  will  introduce  new  modes  with  large 
modal  walues*  Thus*  the  homogenous  etiuationaa  for  instance*  will  have  solu- 
tiona  with  terms  in  the  form  e^  idiere  IM  is  large*  If  X  has  a  large  positlTB 
real  part*  this  term  will  make  the  entire  solution' unstable*  If  X  has  a  large 
negatlwe  real  part*  these  terms  disappesor  rapidly*  the  desired  result*  Further* 
more*  in  this  case  by  proper  machine  procedure  the  effects  of  the  additional 
terms  can  be  made  negligible*  There  is  an  intermediate  case  in  lixieh  rapid 
oseilletiOBs  can  appear* 

These  then  are  the  major  possible  effects  of  the  increase  in  order*  Wo 
refer  to  these  as  the  ^effect*  In  addition*  a  differential  analyser  may 
be  perturbed  so  as  to  chAnge  the  equation  without  increasing  the  order*  These 
we  call  a  errors*  A  differential  analyzer  is  also  dependent  on  stored  infozma** 
tion*  A  perturbation  in  this  stored  information  we  call  a  S  error*  It  is 
possible  to  analyze  the  growth  of  error  in  digital  method  of  solution  in  terms 
of  B  errors*  The  truncation  errors  and  the  round-off  effect  at  each  step 
can  be  considered  as  a  perturbation  on  stored  information*  (Cf*  Rsf* 

The  use  of  digital  check  solutions  requires  a  precise  analysis  of  Ihe  error 
in  the  solution  since  the  solution  must  be  good  enough  to  sexre  as  a  standard 
but  not  so  good  as  to  represent  a  waste  of  oemputation*  Qy  means  of  the  abore 
discussion  of  ^errors*  the  problem  can  be  referred  back  to  ttie  error  at  each 
step*  In  general  this  consists  of  a  truncation  error  and  round  off  errors* 
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let  us  consider  the  truncation  error.  This  will  be  dependent  on  the 
method  of  integration  used.  The  earlier  automatic  sequence  digital  calculators 
had  very  limited  memory  available  and  it  was  urged  that  one  use  methods  like  the 
Runge-Bitta  method  which  are  not  dependent  on  the  storage  of  previously  calcu¬ 
lated  results.  However,  this  argument  does  not  appeal  to  me  exc^t  in  the  case 
of  extremely  limited  storage  and  very  simple  equations  since  the  total  amount 
of  ccqputatlon  is  tripled  in  the  more  general  cases.  Instead  I  would  prefer 
methods  like  the  Ullne  "open"  and  "clos^"  or  simply  "open*  proeeduresi  irtiere 
the  previously  calculated  values  of  the  derivatives  are  utilized  over  again 
to  obtain  higher  order  approximations.  The  additional  storage  required  is 
seldom  significant. 

Among  other  advantages,  one  can  readily  calculate  the  truncation  error 
for  a  p'th  order  open  process  if  the  (p4’2)nd  derivative  is  reasonably  constant 
over  the  interval  Involved  in  the  integration  proeesai  x 


3 


One  can  obtain  a  similar  expression  for  the  truncation  error  for  the  closed 
step, 

B  h*^^  a _ 2. 

9  dt 

80  that  if  one  perfoxns  an  open  step  and  then  a  closed  step,  the  difference 
is  a  fixed  constant  multiple  of  ttie  truncation  error  for  either  step.  One 
can  show  that  there  is  relatively  little  advantage  in  going  beyond  the  first 
closed  step*  (Cf*  Ref.  PJ) 

Since  the  error  for  the  closed  step  is  definitely  aoaller  than  that  for 
the  open  procedure,  this  would  indicate  that  the  prefered  method  might  be  an 
open  and  closed  process  at  each  step.  However,  this  means  that  one  is  doing 
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twio*  the  eork  of  a  single  open  step*  so  that  one  should  ee^pare  the  open  and 
elosed  pvoeedure  constant  with  an  open  procedure  for  one  half  the  intermit 
i.e*  1^  should  be  conpared  with  If  these  quantities  are  coiqparablet  the 

open  nstbed  «lth  the  asaller  interral  is  preferable  since  the  hppotheses  of 
the  discussion  are  aore  strongly  fulfilled* 

I  hare  used  this  procedure  but  with  erery  tenth  step  an  open  followed  by 
a  closed  so  that  one  can  estinate  the  errors*  (Rsf*  [!}) 

Ihe  total  error  can  be  obtained  in  two  ways*  Che  is  the  detaileA  0  error 
■sthod  indicated  above*  On  the  other  hand*  the  S  error  theory  Shows  that  under 
reasonable  general  circumstances  the  total  error  Is  linear  in  the  individual 
truncation  errors*  If  possible*  one  should  therefore  repeat  the  procedure 
with  half  the  interval*  Ibr  if  we  use  half  the  interml*  the  error  is  eut 
doim  to  1/2^  its  fomsr  value  so  that  the  difference  between  the  original 
sslution  and  tts  aolntion  with  half  ttis  interval  is  class  to  the  total  errcr 
for  the  first*  Incidentally*  the  ■extrapolation  to  aero*  proeess  suggested 
by  this  is  based  on  certain  assunptions  concerning  the  order  of  the  error  and 
in  effect  iqproves  the  answer  only  by  an  increase  in  one  in  the  ozder*  flush 
an  inprovcunat  in  the  order  could  be  obtained  directly  by  the  open  integantien 
proeess  in  the  analytic  case  with  less  eo^utation* 

Hoxnal  nunsrical  analysis  is  based  on  the  approziantion  of  arbitrary 
ftinetions  by  polynomials*  Ihis  is  the  reason  for  the  assuaption  aads  above 
that  the  highest  dsrimtive  should  be  essentially  a  constant*  However*  in  the 
differsntinl  equations  which  we  have  considered*  one  can  obtain  an  approaiaate 
esqprsssion  valid  over  a  larger  region  if*  instead  of  using  polynonials*  one 
uses  linear  combinations  of  exponentials*  Fbrthsxmors*  the  *  linearization* 
referred  to  above  seems  to  indicate  that  ezponsntials  would  be  preferable  end 
also  what  siqpensntial  Ametions  to  use* 
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This  of  course  raised  the  problem  of  extending  the  usual  results  of 
numerical  analysis  to  the  situation  in  which  linear  combinations  of  exponentials 
are  used  Instead  of  polynomials.  One  can  regard  the  polynomial  case«  in  fact* 
as  a  special  case  of  the  e3q)onential. 

This  was  basically  the  problem  that  Dr.  Brock  and  I  considered  in  our 
paper  on  the  use  of  exponential  sums  in  integration.  (Cf.  Ref.  £2^)  We 
considered  specifically  the  processes  for  open  and  closed  integration,  emd  by 
the  use  of  certain  detexninants  we  obtained  formulae  which  described  the 
truncation  error  for  various  types  of  integration.  The  method  we  used  for 
evaluating  the  error  is  applicable  to  any  integration  process  which  uses  equally 
spaced  points  and  is  linear  in  the  value  of  the  integrsmd.  Ihese  of  course 
include  such  methods  as  the  open  and  closed  mine  process.  Runge-Bitta.  etc. 
These  formulae  specialize  in  the  case  of  polynomials  to  a  result  Which  is  mors 
convenient  to  use  than  the  usual  derivative  relation* 


62 


Conference  of  Areeml  Nathentleiaae 


BIBlJOQRimY 


[1]  P*  Brock  and  F.  J.  lihi^ayt  "Planning  and  error  analysis  for  the  nuaarieal 
solution  of  a  test  cystem  of  differential  equations  on  the  lai  sequenee  ealeulator" » 
Project  Cyclone  Report  t  ReeTes  Instrument  Corpora tiont  215  >•  91  9t»»  New  York  28 « 

N.  Y.  (1949) 


(21  P.  Brock  and  F.  J.  Uirray«  "The  use  of  exponential  sums  in  step  by  step 
integration* t  Math.  Tables  and  Other  Aids  to  Cooputation*  6  (1952)t  P«  63-78 
and  p.  138-150. 


K.  S.  Miller  and  F.  J.  lAirray.  mathematical  basis  for  an  error  analysis 
of  differential  analyzers*.  Jour,  of  Math,  and  Physics.  ▼.  32  (1953).  P*  136^163. 


[4)  K.  S.  lailer  and  F.  J.  Murray.  "Sxistance  theorems  for  ordinary  differan* 
tial  equations*.  New  York  Uniwersity  Preaa.  N.  Y. .  N.  Y.  (1954). 


(5)  F.  J.  Murray.  "Plaxming  and  error  considerations  for  the  numerical  solu¬ 
tion  of  differential  equations  on  a  sequenee  calculator*.  Mathematical  Tables 
and  Other  Aids  to  Confutation,  r.  4  (1950).  p*  133'’^* 


[61  F.  J.  Murray.  *lrror  analysis  f«r  matbsmatieal  machiaas*.  Trana.  of  the 
N.  Y.  Acad.  Sei..  ▼.  I3  (1951).  P*  168-174* 


^  F.  J.  Murray.  "Mathematical  error  analysis  for  continuous  eonputers*. 
Project  Cyclone  Si^osium  H.  Part  IZ.  Reeres  Inatrummat  Coxpomtiont  0.  8. 
Navy  Special  Oarlees  Center  and  U.  8.  Nawy  Bureau  of  Aaronauties.  (1952) 
p.  139-146.  New  York. 


OEil  J.  Winson.  *The  error  analysis  of  electronic  analogue  eo^mtatian.  linear 
differential  equations  with  constant  coefficients”.  Ph.  D.  Ihesla.  Colviabia 
Unirersity,  (1954)  N.  Y.,  N.  Y. 


